Conductivity of graphene on boron nitride substrates 
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We calculate theoretically the disorder-limited conductivity of monolayer and bilayer graphene 
on hexagonal boron nitride (h-BN) substrates, comparing our theoretical results with the recent 
experimental results. The comparison leads to a direct quantitative estimate of the underlying 
disorder strength for both short-range and long-range disorder in the graphene on h-BN system. 
We find that the good interface quality between graphene and h-BN leads to strongly suppressed 
charged impurity scattering compared with the corresponding Si02 substrate case, thus producing 
very high mobility for the graphene on h-BN system. 

PACS numbers: 72.80.Vp, 81.05.ue, 72.10.-d, 73.22.Pr 



An important recent development in the physics and 
materials science of graphene [H, 0] is the successful fab- 
rication of gated graphene layers on hexagonal boron ni- 
tride (h-BN) substrates 0, 01 • Since h-BN has the same 
hexagonal honeycomb lattice structure as graphene it- 
self with an almost matching lattice constant, the ex- 
pectation has been that graphene on h-BN would have 
much lower disorder than the generic graphene on Si02 
substrates that has almost universally been studied so 
far experimentally. This expectation has, in fact, been 
spectacularly borne out by the recent experiments 0, |3| 
at Columbia University where both monolayer graphene 
(MLG) and bilayer graphene (BLG) on h-BN substrates 
have been shown to have substantially higher (by roughly 
one order of magnitude or more) carrier mobility than 
graphene samples on the standard Si02 substrates [H, 
In fact, the quality of graphene on h-BN, as measured by 
transport experiments, appears to be comparable to that 
of annealed suspended graphene@,[(>], with both systems 
exhibiting clear fractional quantum Hall effects attesting 
to their very high mobility. 

In the current work, we consider theoretically elec- 
tronic transport properties of graphene/BN MLG and 
BLG systems, using the highly successful Boltzmann- 
Kubo-RPA formalism which has earlier been used[l| to 
study graphene transport on Si02 substrates, both for 
MLG [7H9| and BLG [lO| systems, as well as for the 
suspended graphene [ll| system. Our goal is a thor- 
ough quantitative understanding of the specific opera- 
tional features of resistive scattering mechanisms limit- 
ing carrier mobility in graphene on h-BN. By demanding 
quantitative agreement between our calculated graphene 
(on h-BN) transport properties with the corresponding 
experimental data 0, |3| for both MLG and BLG sys- 
tems, we establish the precise role of long-range (e.g. 
charged Coulomb impurities in the environment) ver- 



sus short-range (e.g. point defects, neutral scatterers, 
vacancies) disorder in graphene on h-BN systems. We 
obtain excellent agreement with the experimental data 

13] using very reasonable disorder parameters, estab- 
lishing that the good interface quality between graphene 
and h-BN (e.g. lack of dangling bonds) leads to strongly 
suppressed charged impurity scattering compared with 
the corresponding Si02 substrate situation 0, Q, thus 
providing very high mobility for the graphene on h-BN 
system. The relative suppression of long-range scatter- 
ing compared with the short-range scattering also leads 
to rather nonlinear-looking MLG conductivity as a func- 
tion of gate voltage (i.e. carrier density) for graphene 
on h-BN substrates compared with the Si02 substrates, 
thus explaining the peculiar experimental finding that 
the observed BLG (MLG) conductivity on h-BN sub- 
strates [3, 0] manifests linear (nonlinear) conductivity as 
a function of the gate voltage. Our theory also natu- 
rally explains the weaker observed temperature depen- 
dence of MLG conductivity than the BLG conductivity 
for graphene on h-BN substrate. The actual conductiv- 
ity of MLG/BN or BLG/BN [i,13] is determined by the 
detailed interplay between the long-range and the short- 
range disorder in the relevant system along with the dis- 
tinct screening properties of the graphene carriers as it is 
in the usual graphene on Si02 substrates [H. [tI. Isl. lol 1^ . 

The graphene conductivity a is given by [l| 
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where / — f{sk) is the Fermi distribution function, D{e) 
is the density of states, u/j — dsk/dk is the carrier veloc- 
ity, and T[e) is the transport scattering (or relaxation) 
time which depends explicitly on the effective disorder 
scattering potential 'F' : 
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where Ek is the graphene carrier energy dispersion for 2D 
wave vector k, 'z' is the position of the impurity whose 
concentration is defined by with a denoting the kind 
of impurity (e.g. long-range or short-range), g{9) denotes 
a known chiral matrix element form factor determined by 
the band structure (and is therefore different for MLG 
and BLG), and (k, k') are the incoming and outgoing 
carrier 2D wave vector due to the impurity scattering po- 
tential Vkk'{z). Since the details of the transport theory 
for graphene have been discussed earlier in the literature 
we only make a few comments on the calcula- 
tional aspects of our theoretical results presented in this 
work: (i) the substrate h-BN is characterized by its static 
dielectric constant kbn = 7.0 13|, leading to an effective 
background dielectric constant n = 4.0, which enters into 
the definition of the effective disorder potential; (ii) the 
effective disorder potential 'V entering Eq. (2) is taken 
to be the screened disorder where the screening is by the 
static graphene (MLG or BLG) dielectric function e(g, T) 
which has been calculated earlier in refs. [3,111 respec- 
tively for MLG and BLG; (iii) we include two types of 
disorder in our theory, the long-range disorder character- 
ized by randomly distributed charged impurity centers 
with 2D density of rii located at the graphene-BN in- 
terface and the short-range disorder characterized by an 
effective strength of UdVQ denoting a white-noise delta- 
correlated local disorder. (We emphasize that both long- 
and short-range disorder are necessary for quantitative 
and qualitative understanding of experimental data.) 

The theory is characterized by two parameters, Ui and 
ndVQ, describing long-range and short-range disorder, re- 
spectively. In principle, the effective separation (d) be- 
tween the location of the charged impurity centers and 
the 2D graphene layer could also be an additional physi- 
cally relevant parameter in the transport theory [H, 0, [1] , 
but we put d = throughout this paper, keeping the 
number of free parameters a minimum (only two) and 
assuming that the random charged impurity centers are 
located at the graphene/BN interface as consistent with 
the very high quality of the h-BN crystals used in refs. 
0, 3- I^or obtaining our theoretical transport results 
we have varied the parameters rii and ndV^ arbitrarily 
over a wide range, obtaining the best regression fit to the 
high-density data of refs. [H 0] as shown in Figs. 1 and 
2. (We have also used different values of 'd' using the 
charged impurity separation as a tuning parameters, but 
our results are qualitatively unaffected by an adjustable 
d.) 

We first show our theoretical results valid at "high" 
carrier density (n) defined as n > away from the mini- 
mum conductivity Dirac point regime, where the density 
fluctuations associated with the inhomogeneous puddle 
formation can be safely neglected. In Figs. 1 and 2 we 
show our calculated conductivity (at T = 0), cr(n), as a 
function of the carrier density for a few different values 
of the disorder parameters choosing the parameters such 
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FIG. 1. (Color online) (a) Calculated conductivity of MLG 
(umlg) and BLG {oblg) using the following parameters: 
ni = 12 X 10^" and mVS = 0.7 {eW kf . In (b) [(c)] 
we show the individual MLG [BLG] conductivity limited by 
long-range scattering (oi) and short-range scattering (o"d). o"t 
indicates the total conductivity limited by the both scatter- 
ings. 



that we get essentially exact quantitative agreement away 
from the Dirac point (n > n.i) with the experimental data 
of ref. [i^ for MLG/BN (Fig.l) and BLG/BN (Fig.2) sys- 
tems. In each figure, we present results for both MLG 
and BLG systems for a fixed set of values of the disor- 
der parameters rii (long-range) and UciVq (short-range) 
with the results of Fig. 1 (2) showing quantitative agree- 
ment with the corresponding experimental data for MLG 
(BLG) on h-BN in ref. [3j]. In each figure we present the 
individual conductivity limited by long-range and short- 
range scattering as well as the total conductivity. 

Three qualitative features of our theoretical results in 
Figs. 1 and 2 stand out: (i) for fixed disorder, MLG con- 
ductivity is always larger than BLG conductivity for all 
densities although they approach each other at very high 
density as expected; (ii) the quantitative values of the dis- 
order parameters (i.e. and n^VJ)^) necessary in Figs. 1 
and 2 for obtaining agreement with the experimental data 
for graphene on h-BN substrates are typically much 
(by more than an order of magnitude) smaller than that 
needed for agreement between theory and experiment 
with the corresponding graphene on Si02 substrates (e.g. 
refs. 0,i)— this is particularly true for the charged im- 
purity density n.i which has the remarkably small value 
of 0.3 X 10" - 1.0 X 10" cm-2 for graphene on h-BN 
substrates compared with n.i > 10^^ cm^^ for graphene 
on Si02 substrates (we note that short-range disorder 
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FIG. 2. (Color online) (a) Calculated conductivity of MLG 
(o"mlg) and BLG {(Jblg) using the following parameters: 
m = 4x W^°cm-^ and naVo = 0.2 (eVA)^. In (b) [(c)] 
we show the individual MLG [BLG] conductivity Umited by 
long-range scattering (ai) and short-range scattering {(Jd)- ot 
indicates the total conductivity limited by the both scatter- 
ings. 



characterized by n^Vp^ seems comparable in strength for 
h-BN and Si02 substrates with h-BN having somewhat 
smaller values); (iii) the MLG conductivity results for h- 
BN substrates are much more sublinear than for the cor- 
responding Si02 substrate case clearly establishing the 
much weaker role of long-range charged impurity scat- 
tering in h-BN systems compared with Si02 systems. 

It is easy to show theoretically P, 0, II] using Eq. (2) 
that the charged impurity scattering limited MLG con- 
ductivity af^^^ on h-BN substrates is given approxi- 
mately by crf^^'^ « 25.7(e^//i)(n/ni) whereas the short- 
range scattering limited conductivity is given by cr^^^'^ ~ 
350(e^//i)/(ndVg^) where UdYg is measured in (eVA)^ 
units. Our MLG numerical results for long-range scat- 
tering shown in Figs. 1 and 2 obey these analytical re- 
lations exactly with the net conductivity being given 
by a = {a~^ + cr^ ^)"^- For the BLG on h-BN sub- 
strates, a simple analytic relation can only be derived for 
the short-range scattering-limited conductivity ctq^"^ = 
QQ.l{e'^ /h){n/ndVQ), which is linear in carrier density, 
with 'n' in units of 10"'^^ cm^^ and ridVi^ in units of 
(eVA)^. The long-range disorder leads to a crf^^ ^ n" 
where a w 1 — 1.3 depending on the parameter regime, 
and no simple analytic relationship can be derived ex- 
cept at very low BLG carrier density where the Coulomb 
disorder is effectively completely screened out since the 
BLG screening wave vector becomes much larger than 
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FIG. 3. The temperature dependence of the MLG conduc- 
tivity (a) and the BLG conductivity (b) for several electron 
densities. In (a) [(b)] the parameters of Fig. 1 [Fig. 2] are 
used. 



the Fermi wave vector. For this very low carrier density 
regime (<t; 10^^ cm~^), the charged impurity disorder 
limited BLG conductivity becomes linear in carrier den- 
sity (i.e., a — 1) obeying the approximate relationship: 
(jf^'-^ ^ 15(e^//i)(n/ni). We mention, however, that this 
formula is not useful for n < Ui since density inhomogene- 
ity effects associated with puddles would dominate close 
to the charge neutrality point. 

The agreement, using reasonable values of disorder pa- 
rameters, between theoretical results presented in Figs. 1 
and 2 with the experimental data 0, 3| of the Columbia 
group indicates that graphene on h-BN indeed has sub- 
stantially lower long-range Coulomb disorder in its envi- 
ronment than graphene on Si02 substrates, most likely 
due to the high-quality graphene/BN interface without 
any dangling bonds as already speculated in ref. A 
direct consequence of this reduced Coulomb scattering is 
the manifestly sublinear a{n) observed in the MLG/BN 
system to be contrasted with the linear cr(n) in the 
MLG/Si02 system P, Q except at very high densities. 
We note that our theory indicates a direct way of estimat- 
ing the strength of both long- and short-range disorder 
from the high-density MLG/BN a{n) data by obtaining 
the slope da/dn at high-density (which gives Ui) and by 
obtaining the intercept of the high-density cr{n) extrap- 
olated to n — > 0, which gives UdVQ. 

In Fig. 3, we present our theoretical results for the 
temperature dependence of the MLG (Fig. 3) and BLG 
(Fig. 4) conductivity, a{n, T), on h-BN substrates. These 
results are again valid (similar to those in Figs. 1 and 
2) at high carrier density (n > n^) where density in- 
homogeneity effects are negligible. All phonon effects 
Il7| are neglected here with the temperature depen- 
dence arising entirely from the temperature dependence 
of the screening function and the energy-averaging asso- 
ciated with the finite-temperature smearing of the Fermi 
surface [9|. The first effect ("screening") produces weak 
metallic temperature dependence (i.e. a decreasing with 
increasing T) since screening weakens at higher temper- 
atures whereas the second effect ("thermal averaging") 
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FIG. 4. (Color online) Calculated conductivity using effec- 
tive medium theory (solid lines) for both MLG (umlg) and 
(o"slg)- Dashed lines represent the transport results calcu- 
lated using Boltzmann transport theory. In (a) and (b) we 
use the parameters corresponding to the Fig. 1 and 2, re- 
spectively. Note that the effective medium theory becomes 
important only at low (n < rii) carrier density where inhomo- 
geneous puddle formation becomes significant. 



produces weak insulating temperature dependence (i.e. a 
decreasing with increasing T) . Although the temperature 
dependence is weak, as is obvious from Fig. 3, the the- 
oretical behavior of a{T) is qualitatively consistent with 
the experimental observations Q away from the Dirac 
point: (i) MLG manifests weak metallic T-dependence, 
and (ii) BLG manifests weak insulating T-dependence. 
Experimentally, both systems manifest insulating a{T) 
at the Dirac point where density inhomogeneity effects 
dominate but at higher density in> rii) our results 
are consistent with experimental finding of ref . Q . 

Finally, we consider in Fig. 4 the low-density trans- 
port in graphene/BN systems, where the results (valid 
for n > rii) shown in Figs. 1-3 do not apply. At low car- 

for 
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rier density n (< n^), which is very low (~ 10""cm" 
the graphene/BN system because of its extremely weak 
Coulomb disorder, the graphene layer is known [l|,[l,[3- 
22| to break up into inhomogeneous puddles due to the 
failure of screening, and thus the naive Boltzmann-Kubo- 
RPA transport theory, which explicitly assumes a homo- 
geneous carrier density, is no longer valid since the den- 
sity fluctuations become strong. To demonstrate the ef- 
fect of puddle formation on graphene/BN transport prop- 
erties, we have carrier out an effective medium theory 
lOl |23| calculation of transport using rirms = ni, where 
rirms is the root-mean-square fluctuation in the carrier 
density due to the puddles induced by the charged impu- 
rities. Microscopic self-consistent calculations l3, 3 23 1 
show that rirms ~ rii is a reasonable qualitative approx- 
imation for the density inhomogeneity around the Dirac 
point. Our T = effective medium theory transport re- 
sults (using the Boltzmann-Kubo-RPA transport formal- 
ism) are shown in Fig. 5 for both MLG/BN and BLG/BN 
systems. The most important features of Fig. 5 are: (i) 
the high density {n > rii) results shown in Figs. 1 and 
2 remain valid; (ii) near the Dirac point (for n < rii), 
(j{n) saturates with a non-universal disorder-dependent 



minimum conductivity [l|, H, whose value is roughly 
given by 2 — 10 (e^ /h) consistent with the experimental 
observations 

We have also carried out our effective medium theory 
calculation at finite temperature to include the puddle 
effects on a{n,T). These results (not shown) agree with 
the results shown in Figs. 3 and 4 at high densities {n > 
rii), but for low carrier densities (n < rii) we qualitatively 
recover the experimentally observed strongly insulating 
cr(r) induced by the puddles. 

This work is funded by US-ONR. 



[1 
[2 

[3: 

[5: 
[6: 

[r 

[8 
[9 

[lo; 

[11 

[12 
[13 

[14 
[15 
[16 
[17 



S. Das Sarma, S. Adam, E. H. Hwang, and E. Rossi, 
.arXiv: 1003.4731 (Rev. Mod. Phys. in press). 
N. M. R. Peres, Rev. Mod. Phys.82, 2673 (2010); A. H. 
Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, 
and A. K. Geim, Rev. Mod. Phys. 81, 109 (2009). 
C.R. Dean, A.F. Young, I. Meric, C. Lee, L. Wang, S. 
Sorgenfrei, K. Watanabe, T. Taniguchi, P. Kim, K.L. 
Shepard, and J. Hone, Nature Nanotech. 5, 722 (2010). 

C. R. Dean, A.F. Young, P. Cadden-Zimansky, L. Wang, 
H. Ren, K. Watanabe, T. Taniguchi, P. Kim, J. Hone, 
and K.L. Shepard, arXiv: 1010. 1179 

F. Ghahari, Y. Zhao, P. Cadden-Zimansky, K. Bolotin, 
and P. Kim, arXiv:1010.1187 

W. Bao, Z. Zhao, H. Zhang, G. Liu, P. Kratz, L. Jing, 
J. Velasco, D. Smirnov, and C. N. Lau, Phys. Rev. Lett. 
105, 246601 (2010). 

E. H. Hwang et al., S. Adam, and S. Das Sarma, Phys. 
Rev. Lett. 98, 186806 (2007); Y.-W. Tan et al, Y. Zhang, 
K. Bolotin, Y. Zhao, S. Adam, E.H. Hwang, S. Das 
Sarma, H. L. Stormer, and P. Kim, Phys. Rev. Lett. 99, 
246803 (2007). 

S. Adam, E. H. Hwang, V. M. Galitski, and S. Das 
Sarma, Proc. Natl. Acad. Sci. USA 104, 18392 (2007); 
J. H. Chen et al, C. Jang, S. Adam, M. S. Fuhrer, E. D. 
Williams, and M. Ishigami, Nature Phys. 4, 377 (2008). 
E. H. Hwang and S. Das Sarma, Phys. Rev. B79, 165404 

(2009) . 

S. Das Sarma, E. H. Hwang, and E. Rossi, Phys. Rev. 
B81, 161407 (2010). 

K. Bolotin, K. Sikes, Z. Jiang, G. Fudenberg, J. Hone, 
P. Kim, and H. Stormer, Solid State Commun. 146, 351 
(2008); S. Adam and S. Das Sarma, Solid State Commun. 
146, 356 (2008). 

S. Xiao, J. H. Chen, S. Adam, E. D. Williams, and M. S. 
Fuhrer, Phys. Rev. B 82, 041406 (2010). 
R. Geick, C. H. Perry, and G. Rupprecht, Phys. Rev. 146, 
543 (1966); P. J. Gielisse, S. S. Mitra, J. N. Plendl, R. 

D. Griffis, L. C. Mansur, R. Marshall, and E. A. Pascoe 
Phys. Rev. 155, 1039 (1967). 

E. H. Hwang and S. Das Sarma, Phys. Rev. B 75, 205418 

(2007) . 

E. H. Hwang and S. Das Sarma, Phys. Rev. Lett. 101, 
156802 (2008). 

E. H. Hwang and S. Das sarma, Phys. Rev. B77, 115449 

(2008) . 

D. K. Efetov and P. Kim Phys. Rev. Lett. 105, 256805 

(2010) . 



5 



[18] E. H. Hwang and S. Das Sarma, Phys. Rev. B82, 081409 

(2010). 

[19] E. Rossi and S. Das Sarma, Phys. Rev. Lett. 101, 166803 
(2008). 

[20] J. Martin, N. Akerman, G. Ulbricht, T. Lohmann, J. H. 
Smet, K. von Klitzing, and A. Yacobi, Nat. Phys. 4, 144 



(2008). 

[21] Y. Zhang, V. W. Brar, C. Girit, A. Zcttla, and M. F. 

Crommic, Nature Phys. 5, 722 (2009). 
[22] A. Deshpande, W. Bao, Z. Zhao, C. N. Lau, and B. J. 

LeRoy, App. Phys. Lett. 95, 243502 (2009). 
[23] E. Rossi, S. Adam, and S. D. Sarma, Phys. Rev. B 79, 

245423 (2009). 



